# Plot frequency response
import control    #pip install control
import matplotlib.pyplot as plt
import math
import numpy as np
sysTfList=[]
sysTfList.append(control.tf([1],[1,0.4,0]))
mag,phase,omega=control.bode_plot(sysTfList,omega=np.logspace(start=-2,stop=1,num=200),dB=False,deg=True,plot=True)
plt.show()